Neutrino oscillations in a model with a source and detector 
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We study the oscillations of neutrinos in a model in which the neutrino is 
coupled to a localized, idealized source and detector. By varying the spatial 
and temporal resolution of the source and detector we are able to model the 
Q_i' full range of source and detector types ranging from coherent to incoherent. 

O ' We find that this approach is useful in understanding the interface between 

the Quantum Mechanical nature of neutrino oscillations on the one hand and 
the production and detection systems on the other hand. This method can 
easily be extended to study the oscillations of other particles such as the 
neutral K and B mesons. We find that this approach gives a reliable way to 
treat the various ambiguities which arise when one examines the oscillations 
from a wave packet point of view. We demonstrate that the conventional 
oscillation formula is correct in the relativistic limit and that several recent 
claims of an extra factor of 2 in the oscillation length are incorrect. We also 
demonstrate explicitly that the oscillations of neutrinos which have separated 
spatially may be "revived" by a long coherent measurement. 
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I. INTRODUCTION 



The flavour oscillations of particles are a fascinating demonstration of quantum mechan- 
ics in the macroscopic world. Flavour oscillations can generically occur when the states 
which are produced and detected in a given experiment are superpositions of two or more 
eigenstates which have different masses. The oscillations of K and B mesons have been 
observed experimentally Q and have been used to place stringent constraints on physics 
beyond the Standard Model. If neutrinos are massive, they too may oscillate, and this could 
lead to the resolution of the well-known solar neutrino problem More recently, the 

discussion of particle oscillations has been extended to include supersymmetric particles in 
super symmetric extensions of the Standard Model ||. 

The phenomenon of particle oscillations has been studied extensively and is generally 
thought to be very well understood. There nevertheless remain several subtle issues which 
continue to cause some confusion. The key to a complete understanding of any such issue 
lays in treating correctly the necessary interplay between the "classical" and "quantum" 
natures of the particles which are interfering to produce the oscillations. Thus, for example, 
the interference effect itself is purely "quantum" in nature (it requires that the particles 
be described by waves), and yet the resulting oscillations in space are only observable if 
the particles are sufficiently localized in space ||. This example highlights the fact that 
any discussion of particle oscillations implicitly assumes that the mass eigenstates which are 
interfering to produce the oscillations are described by some sort of wave packets. 

Despite the success of the wave packet approach in clarifying many aspects of the phe- 
nomenon of particle oscillations |6|-|9[], the approach is not without its own difficulties. The 
results of a given calculation will depend, for example, on the details of the initial mass 
eigenstate wave packets (including their shape, spectrum and relative normalization). One 
particularly difficult problem which arises is the conversion of the final time-evolved wave 
packets into an experimentally observable quantity: since it is generally the flux of particles 
which is measured in an experiment, one is required to calculate a current density rather 
than a probability density^ The difference between a current density and a probability den- 
sity, at least naively, involves a factor of the velocity v, which is very significant if the mass 
eigenstates have quite different masses. Thus, if one would calculate the probability density 
at the detector and integrate it over time, the resulting expression would have factors of 1/v 
pre-multiplying the various terms, leading to an enhancement of the terms corresponding to 
heavier mass eigenstates. In the case of neutrinos, as was noted in [§-10], this would skew 



the usual oscillation formula quite dramatically if one of the mass eigenstate neutrinos was 
non-relativistic. Efforts to construct an appropriate current density which retains the nec- 
essary wave packet features have had mixed success. A calculation in the kaon case appears 



to give reliable results ||11|| , but it can be shown that unphysical effects arise if one attempts 
to define a suitable current when the mass eigenstates have very different masses JTj]. 

There is another very striking apparent "ambiguity" which arises if one does not treat 



Wave packet calculations lead naturally to expressions for the probability density, which are 
appropriately integrated over space, not time. For oscillations in space one wants a quantity which 
is appropriately integrated over time, i.e., a current. 
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the delicate interplay between the classical and quantum natures of the particles correctly. 
In this case the "ambiguity" leads to an alleged error by a factor of "2" in the calculation 
of the oscillation length In order to understand the source of the ambiguity, we follow 



the discussion given by Lipkin . Suppose we consider the oscillations in time of a system 



for which the initial and final states are not eigenstates of the free Hamiltonian. The phase 
of the interference term will then be given by 4>(t) = (Ei — E 2 )t, where Ei=(p 2 + m 2 ) 1 ^ 2 . 
Detectors do not measure oscillations directly as a function of time, however, so one needs 
somehow to convert this expression into an oscillation in terms of space. We may then, as is 
conventionally done, set x=vt, with v=2p/ (E 1 + E 2 ) representing a sort of average velocity. 
We then obtain the following phase in terms of x, 

= K : mi)I . (i.i) 

2p 

This is the conventional (and correct) result for the phase difference. Let us now attempt to 
incorporate a classical aspect of the problem and argue that since the two mass eigenstates 
travel at different speeds, they will arrive at the detector at differenttim.es t\ and t 2 , related 
by x=pti/ Ei=pt 2 / E 2 . Taking the phase of the interference term to be 4>(ti,t 2 )=E 1 ti — E 2 t 2 , 
we then obtain 

*,„(*) = W -"■*>* , (1.2) 
p 

which differs from the conventional phase difference, Eq. ( |1.1|) , by a factor of two. This 
result, were it correct, would indeed be rather remarkable. 

The first resolution of this ambiguity was given by Lipkin [|T4]] (see also Ref. [15[]) who 
argued, on physical grounds, that the energies (rather than the momenta) of the two mass 
eigenstates should be set equal. In this case the oscillations are described in terms of 
distances directly (since it is the momenta of the two eigenstates which differ) and the 



"correct" oscillation formula is obtained. Lowe et al |16| and Kayser |I7j have extended 
this discussion and have argued that the key to avoiding ambiguities is to ensure that one 
evaluates the wave functions of the mass eigenstates at precisely the same space-time point. 
That is, even though classically the mass eigenstates will arrive at the detector at different 
times, quantum mechanically the wave functions corresponding to different space-time points 
cannot interfere. Indeed the analysis leading to the expression in Eq. ( |1.2| ) involves interfering 
the wave functions for the two mass eigenstates at the same position but different times and 
hence gives the incorrect result. This issue is, in fact, quite subtle. For example a long 
coherent measurement in time may be used to "revive" particle oscillations even after the 



mass eigenstate wave packets have completely separated spatially [|IBJ.0 Thus, wave packets 
arriving at the "classically" separated times t\ and t 2 - and having negligible overlap in their 



2 This behaviour is analogous to what happens when a high Q oscillator gets hit by two successive 
pulses. The first pulse sets the oscillator in motion, causing it to oscillate for a time determined 
by its Q value. If the oscillator is still oscillating when the second pulse arrives, the resulting 
oscillations will exhibit an interference pattern. 
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wave packets - may still interfere to give rise to oscillations. There is then some sense in 
which wave functions corresponding to different space-time points may interfere. 

In light of the issues presented above, it is our view that a proper treatment of the 
Quantum - Classical interface of particle oscillations should incorportate the source and 
the detector as key components of the system. In this paper we present a simple model 
for a particle source-detector system which addresses many of the above issues in a very 
natural and self-consistent way. We shall, for concreteness, consider the case of neutrino 
oscillations, but our approach could easily be adapted to other situations. The source and 
detector will be modeled by simple harmonic oscillators which are de-excited or excited by 
emitting or absorbing neutrinos of a given flavour. (Two-level "fermionic" source - detector 
systems could also be considered.) Having defined the model, it will be straightforward to 
calculate the oscillation probability as a function of the distance between the source and 
detector. The resulting expressions will be found to exhibit all of the known "wave packet" 
characteristics in the relativistic limit, but will also give useful insight into cases in which one 
or more of the mass eigenstates is non-relativistic. In particular, we will find no evidence for 
the enhancement of non-relativistic neutrinos which can occur in conventional wave packet 
analyses. Including the source explicitly in the calculation gives the added benefit that the 
characteristics of the initial wave packets corresponding to the various mass eigenstates are 
completely determined by the characteristics of the source and need not be put in by hand. 



Our approach is similar in spirit to the calculations in Refs. ||T9| , pO| , but is more transparent 
due to the simplified model which we consider. (See also Ref. |21| for a similar calculation 
performed within the context of elementary quantum mechanics.) One advantage of our 
simplified approach is that the dependence on the time-resolution of the detector is very 
clear. This will allow us to verify explicitly that a long coherent measurement in time may 
be used to revive the oscillations of particles whose wave packets have separated spatially^ 
We shall also settle the issue of the "factor of 2" (hopefully) once and for all. 

We begin in the next section by analyzing a simple model in which the neutrino is 
described by a complex scalar field. This field is coupled to two localized simple harmonic 
oscillators, representing the source and detector. Modeling the neutrino by a complex scalar 
field allows for a simpler and more complete evaluation of physical quantities than if a spinor 
field is used. In Sec. [II A| we study the case of a single neutrino species coupled to the source 
and detector. This allows for a careful analysis of the efficiency of our system at producing 
and detecting neutrinos of different masses. Although no oscillations are possible in this 
case, this calculation will be essential in interpreting the results when neutrino oscillations 



are present. In Sec. [II B| we couple several neutrino fields to the source and detector. This 
gives rise in a natural way to oscillations (as a function of the distance between the source 
and detector) in the probability for the source to decay and the detector to be excited. These 
are of course "neutrino oscillations." Sec. [II Q contains a brief analysis of the non-relativistic 
case. We then extend our analysis in Sec. |TTT| to a more realistic model in which the neutrinos 
are described by Dirac spinor fields. These results are compared to the ones with a complex 
scalar field. We conclude in Sec. [TV] with a summary and discussion of our results. 



3 A recent paper has also demonstrated this effect explicitly 22 . 
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II. A MODEL FOR A NEUTRINO SOURCE AND DETECTOR 



The idea of using an idealized detector to clarify physically measurable quantities in 
Quantum Field Theory has been used extensively in the analysis of Quantum Fields in 
non-inertial frames and in gravitational backgrounds [p3 |. In our idealized model, we have 
chosen to couple the neutrino field to two harmonic oscillators, one representing a neutrino 
"source," and the other representing a neutrino "detector." The neutrinos are first taken to 
be complex scalar fields which simplifies the calculations considerably.^ 

The physical picture which we have in mind is the following: we imagine our "source" and 
"detector" to be microscopic on the scale of some macroscopic "bulk" source and detector, 
but to also be very massive compared to the energy of the exchanged neutrino (so that 
the dynamical degrees of freedom of the source and detector may be ignored). Thus, for 
example, the source (detector) could represent some nucleus inside a bulk sample which 
undergoes beta decay (inverse beta decay). The spatial "widths" of the source and detector 
in our calculation are then widths appropriate to, say, nuclear or atomic dimensions. In 
principle, the oscillation probability which we calculate here should subsequently be averaged 
incoherently over the physical dimensions of the macroscopic source and detector, although 
we do not perform this average. If the size of the macroscopic source and detector are 
much smaller than the neutrino oscillation length (which they need to be in order to observe 
oscillations), then this averaging would have only a small effect. 

The interactions at the source and detector will be made explicitly time dependent so 
that they may be turned "on" and "off." This is in keeping with our physical picture. In 
general a real (microscopic) source or detector will be in an environment which is "noisy," 
so that the coherent emission or absorption of a neutrino gets cut off after some time due 
to the interactions of the source or detector with its surrounding environment |9,|l8|. The 



amount of time which the model source or detector spend being "on" is then related to the 
"coherence time" of the physical source or detector.0 

Our calculation proceeds as follows. We first write a Lagrangian which couples the source 
and detector to the neutrino field. In the initial state of the system, the source is in its first 
excited state (ready to emit a neutrino) and the detector is in its ground state. We then 
calculate the probability that at some time far in the future the source is found to be in its 
ground state and the detector in its first excited state. The model will be constructed in 
such a way that this interaction will correspond to exactly one neutrino being exchanged 
between the source and detector (to first non- vanishing order in perturbation theory.) In 
this approach, then, the neutrinos themselves are not observed, but are simply the exchange 
particles in the source-detector interaction. 



The main drawback of this approach is that it ignores the neutrino's spin and the characteristic 
V — A nature of neutrino interactions. 

5 The explicit turning on and off of the source and detector violates energy conservation microscop- 
ically, but that is natural since the interactions of the source and detector with their respective 
environments involve the exchange of energy. If we choose to look at the source or detector in 
isolation, this exchange of energy appears as energy non-conservation. 
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A. A Single Species of Neutrino 

To describe our model, we begin with a single complex scalar field <f>{x) and two oscillators 
qi(t) and q2(t) describing the source and detector, respectively. The action for our model is 
given by 



S = Jd 4 x (£° + £ int ) + J dtL° q , 



where 



-0'(X 



□ + m 2 ) cf)(x 



(2.1) 

(2.2) 
(2.3) 

(2.4) 

The functions €i(t) are explicit functions of time which allow us to "turn on" and "turn off" 
the interactions, and the functions fta(x) (/i 2 (x)) are smooth functions of x which vanish 
outside the source (detector). 

We quantize the free fields in the usual way, requiring 



C 



int 



q{(t) gi (t) - n{qt(t) qi (t) + ql(t)q 2 (t) - fi&J (*)&(*), 
-e x {t) ((P\x)q 1 (t)h 1 ( X ) + <j>(x)ql(t)hZ(x) 
-e 2 (t) U\x)q 2 (t)h 2 (x) + <P(x)ql(t)h* 2 (x) 



(x,t),vr(y,t)] =^ 3 (x-y) 
h(t),Pi(t)] =i. 



(2.5) 
(2.6) 



All other commutators are taken to vanish. The field operators may then be expressed in 
terms of creation and annihilation operators as follows 



where 



dk (a(k) 
2^( A 

dk = 



( '— + b ] (k)e- 

iO.it _|_ ^tgiHit ' 



d 3 k 



(2tt) 3 2£ 

and where the annihilation and creation operators satisfy the commutation relations 
a{k),a\k')\ = \b{k),tf{k')] = (2vr) 3 2£5 3 (k - k'), 



BuB] 



2a 



(2.7) 
(2.8) 

(2.9) 

(2.10) 
(2.11) 



We interpret a){k) and a(k) in the usual way as the operators which create and annihilate, 
respectively, a neutrino state with four-momentum k. W{k) and b(k) act similarly with 
respect to the anti-neutrino states. The operators A\ and A; and B\ and Bi interpolate 
between the energy levels of the harmonic oscillators.^ 



6 Note that we have allowed the to be complex. Had we not done this, the source and detector 
would have exchanged both neutrinos and anti-neutrinos. 
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We take as our initial state 



-oo) = |O;l;O) = |O) ®|l)i®|O) 5 



(2-12) 



in which 



|i)i = 4|o)* 



(2.13) 



represents the first excited state of the oscillator i and in which |0)^ is the neutrino vacuum 
state. We wish to calculate the amplitude for the process in which the source de-excites to 
its ground state and the detector is excited to its first excited state. That is, 



A = (0; 0; l\s, oo) = (0; 0; l|Texp 



—i 



H s (t')dt' 



-oo), 



(2.14) 



in which H s represents the Hamiltonian in the Schrodinger picture. The modulus squared 
of this amplitude is the probability for the transition to take place. 

We shall assume the couplings in the interaction Hamiltonian to be sufficiently small that 
the amplitude in Eq. (|2.14|) is always much less than unity. This is of course always the case 
in the real-world situation which we are attempting to model - neutrino interactions are so 
weak that perturbation theory is always valid. It is then straightforward to evaluate (|2.14|) 
using standard techniques to obtain to leading order and up to an over-all unobservable 
phase, 



A 



-^(0;0;1\T 



tfdt»H*(f)H*(f) 



|0;1;0), 



(2.15) 



where H? t (t) refers to the interaction Hamiltonian evaluated in terms of the free fields in 
the Heisenberg picture at time t. The above expression may be evaluated explicitly in terms 



of neutrino propagators (|P9| , pO! for arbitrary turn-on/off functions ti{t). We find it simpler, 
however, to "design" the turn on/off functions so that the source and detector are never on 
at the same time, and, furthermore, so that the source always turns on first and only then the 
detector. (This avoids the unphysical situation in which the detector emits an anti-neutrino 
which is subsequently absorbed by the source. The amplitude for this process would in 
any case be very small since it violates energy conservation.) Under this assumption only 
one of the time-orderings in the propagator gets picked up and A may be evaluated using 
Eqs. 0, (pD, flnOD , ( gT[P and ( gTg ) to obtain 



A = -(0;0; 1| / dt'dt"d 3 x'd 3 x"e 1 {t')e 2 {t") 

= - j dt'dt"d 3 x'd 3 x"dke 1 {t')e 2 {t")h 1 (x.')h* 2 (x.") 

x exp [-i(E - n 2 )t" + i(E - n x )t' + ik ■ (x" - x')] 



(2.16) 



(2-17) 



Since the amplitude is proportional to (0\4>(x")4>^ (x')\0) , it is clear from Eq. (|2.7| ) that this 
interaction corresponds to the creation and subsequent annihilation of a single neutrino. 

In order to proceed further we choose hi, h 2 and e% to be Gaussians since this allows 
many of the integrals to be evaluated exactly. Setting 
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h 1 ( X ) = (V2^* Xi y 3 e^ 2/2 <, 

/ i2 (x) = (v / 2^.r, 2 )" 3 e- | ™ |2/2 <, 
e 1 (t) = e° 1 e~ t2/2 < 

we obtain 

i{E-n 2 )t" -^{E-n x fal 

(2.21) 
where 

k = VE 2 - m 2 . (2.22) 

Before choosing an explicit form for e^it)^ which determines characteristics of the detec- 
tor, let us make a few observations regarding the above expression for the amplitude. First 
of all, for large Xp, the amplitude decreases like x^ 1 so that the probability falls like x^ 2 , as 
expected on geometrical grounds in three dimensions. At the origin, however, the amplitude 
does not diverge (despite the 1/xd factor), due to the sine function in the integrand. A sec- 
ond observation is that conservation of energy at the source and of momentum at both the 
source and detector are governed by the relative sizes of a tl , u Xl and a X2 . This situation is 
in accordance with the uncertainty principle (and is in fact necessary, as discussed above, in 
order to observe oscillations). In general, neither energy nor momentum need be conserved 
exactly if the source and detector are localized in space and time. The specific set-up which 
we have chosen favours energies close to the energy of the excited source, Q\, and momenta 
close to zero. This latter point is due to the fact that our souce and detector have no dy- 
namical degrees of freedom - they cannot recoil when a neutrino is emitted or absorbed - 
and thus the neutrino gets all of its momentum from the uncertainties in the positions of 
the source and detector. In order to avoid the problem that low momenta are favoured, we 
shall typically choose to set a tl ^ > cr Xl 2 m our numerical work below.Q When several neutrino 
fields are coupled to the source and detector, this will mean that the energies of the mass 
eigenstates will be approximately equal, while their momenta will be determined by their 
energies. Furthermore, the sizes of the neutrino wave packets will then be determined more 
by the amount of time for which the source emits an uninterrupted wave-train than by the 
localization of the source-field interaction in configuration space. In Sec. PH , when we extend 
our analysis to fermionic neutrinos, we will allow the source to decay by emitting both a 
neutrino and its associated lepton. In this case the neutrino's momentum will no longer be 
centered about k^O. 



(2.18) 
(2.19) 
(2.20) 



A = - 



2?reK 
v Ah 2 xd 



dt"e 2 (t") 
sin(kxi)) 



dE exp 



7 This is a "trick" which we use to get sensible results, but it is also not unreasonable on physical 
grounds. According to the discussion in Ref. [18|, for example, this condition is satisfied by several 
orders of magnitude if a x is taken to be on the order of nuclear sizes. The reader is also referred 



to the discussion of Lipkin [14], where this same point is emphasized. 



Let us now study the system as a function of the coherence time of the detector. At one 
extreme we can imagine that a given (microscopic) detector is turned on for the entire time 
that the neutrino "wave packet" passes by. This is an ultimately "coherent" detection event. 
Another possibility is that a given microscopic detector turns on and off without sampling 
the entire wave packet. In order to model the former scenario we use a simple step function 
for 62, while for the latter case we use a gaussian: 

^W^efyfa-tMt-h), (2.23) 
ef uss (0 = e° 2 e- {t - tD)2/2 <. (2.24) 

The step function detector turns on abruptly at time t\ and off again abruptly at time t 2 
(with t\ and t 2 chosen such that the entire wave packet passes by while the detector is on), 
while the gaussian detector turns on and off gradually at a time centered around to- Since 
the coherent (step function) detector "catches" the entire wave packet, there is no need to 
integrate the resulting expression for the probability over time. This is not the case for the 
incoherent (gaussian) detector: since each individual microscopic detector sees only a piece 
of the wave packet, one must sum incoherently over all of the microscopic detectors in order 
to correctly model the response of the bulk detector. If there are many microscopic detectors 
in the bulk detector, this incoherent sum corresponds to integrating the expression for the 
probability over to- 
ll is straightforward to evaluate the amplitudes for both types of detectors and we obtain 



= m - tl) j m iE ^JS_n exp 



- l -{E-Q 2 ){t 1 +t 2 ) 



'-\{E-^fal-\{E-Q 2 fal 



sin(fcr D ), (2.25) 



-4 g auss = Ny/2na t2 f dE exp 

Jm 



1 

2 

where 



~k 2 (al+al 2 )-i(E-n 2 )t D 



sm(kx D ), (2.26) 



N = 1 2 fl ( 2 27) 

It is not possible in general to obtain analytic closed-form solutions of these integrals, but 
they are simple to evaluate numerically. In so doing, we obtain exact (to second order in 
perturbation theory) solutions to the problem which we are studying, including all effects 
due to the spreading of the neutrino wave packets. Alternatively, we may, in some cases, find 
reliable approximations for these integrals. Such is the case for the step function detector 
if the neutrino's mass is not too close to the production and detection thresholds (m<^f2j — 
1/otJ. Taking t\ to be a time before the first bit of neutrino flux arrives at the detector and 
t 2 to be a time after the entire neutrino wave packet has passed (formally, t 2 ^oo), we find 



lim A st ep(xD,ti,t 2 ) ~ -iNnexp 

<2— >0O 



ikxo - ^(^2 - ^i) 2 of x 



~(nl-m 2 )« + 



(2.28) 
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where k=(Q 2 — m 2 ) 1 / 2 . (The details of this calculation may be found in Appendix 0.) Note 
that the coherent detector "picks out" momenta corresponding to the energy 
The modified "probability" associated with the above amplitude is given by 

V stcp (x D ) = \A stcp (x D )\ 2 /N\ (2.29) 

where we have dropped the implicit dependence on t\ and t%. We have also divided through 
by N 2 because the value of that constant (including the fall-off not really of interest 

to us since in any calculation of the oscillation probability N 2 always factors out. Taking 
the ratio of this probability for two different values of the mass reveals that the system is 
more efficient at producing and detecting higher-mass neutrinos: 



P ste p(m;x D ) 

^ — 7 — ~i V - ex P 

V s t e p{m';x D ) 



(m 2 -m' 2 )« + a, 2 2 ) . (2.30) 



The mass-dependence of the source/detector system arises due to the fact that our source 
and detector favour neutrino states with momenta close to zero. This feature was predicted 
already in the discussion following Eq. fl2.21|) and is due to the fact that the source and 



detector in our model cannot "recoil" and thus the neutrino gets all of its momentum from 
the uncertainty in the positions of the source and detector. Thus the upper limit on the 
neutrino's momentum is given by A; max ~l / a xi 2 . Note that this preference for non-relativistic 
neutrinos is essentially a quirk of our model and should not be viewed as a physical effect. 
The mass-dependence of the system can be minimized by setting a xi 2 to be much less than 
£li2- 111 suc h cases, the step function detector becomes nearly "ideal;" that is, it detects 
neutrinos of different masses with nearly the same efficiency. 

We now turn to the gaussian detector and define a modified probability in analogy with 
Eq. ( gg§ 



V ganss (x D ,t D ,a t2 ) = \A ganss (x D ,t D ,a t2 )\ 2 /N 2 . ( 2 -31) 

This expression gives the probability that a given microscopic detector - turned on for a 
time <7t 2 centered around the time to - is excited. We need to convert this expression into 
one giving the probability that the bulk detector "detects" the neutrino (i.e., that one of 
the micropscopic detectors is excited.) We assume that the bulk detector is "on" for all 
t£>>0 - in the sense that at any given time many of the microscopic detectors are "on" - 
but that the microscopic detectors themselves turn on and off randomly, so that the number 
which are "on" at any given time is roughly constant. Then the probability that the bulk 
detector "detects" the neutrino is proportional to the integral of Eq. (|2.31|) over £d.0 We thus 



Consider first a simpler case in which there are N detectors, turning on and off at times centered 
about ti<t2<- ■ .<t]y. Each of them has probability e to detect the neutrino, but only if one of the 
previous detectors has not already detected it. Then the probability that none of them detects the 
neutrino is (1 — e) , that the last one detects it is (1 — e) N ~ 1 e, that the second last one detects it is 
(1 — e) N ~ 2 e, and so on. The probabilities for the N+l distinct possibilities sum to unity, as required. 
The probability that the neutrino is detected is then 1 - (1 - e) N =Ne - N\e 2 /(N - 2)!2! + . . ~iVe 
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refer to this type of bulk detector as an "incoherent" detector, since we sum the probability 
incoherently over different times. 

The time integral of Eq. ( |2.31| ) may actually be done explicitly. Let us define the following 
(unnormalized) time-integrated probability 



/•oo 

"^mcoh^D? 0"t 2 ) = / dt£)V gaiUSS (x£), to, Ct 2 
J 



(2.32) 



Since the integrand is symmetric under tp— >— tp, we may formally extend the integration 
to negative infinity and divide by two. The time integral then reduces to a delta function 
in energy and allows us to perform one of the energy integrals. As a result, we obtain 



V incoh (x D , a t2 ) = 2n 2 al 2 I dEexp —{E - Vt^al - (E - fi 2 ) V 
-k 2 {a 2 Xi +v 2 X2 )} sm 2 {kx D ). 



(2.33) 



If m^Cfij — l/<7 ti and cr tl 2 ^>a xi 2 (the latter condition is always assumed) then we may 
approximate the above expression by setting sin 2 (fc££))~l/2 to yield 



^mcoh (•££> ; 



5/2 2 
*2 



7T ' a. 



1/2 



cxp 



2/2 i 2 \ 



x exp 



a 



ti + 4 + a lr + °i 



^2.34) 



Thus, under the above conditions the "incoherent" gaussian detector has the same mass- 
dependence as the step function detector does (c.f. Eq. Q2.30| )): 



Vmcoh(m;x D ,a t2 
V incoh {m'; x D , a t2 



exp 



K-^)« + <) 



^2.35) 



This fact is rather remarkable and shows again that it is correct to perform the time integral 
in Eq. (Q). 

Fig. [I] shows a plot of the time- integrated probability, Eq. fl2.33p , as a function of the 
neutrino mass for two different values of u xi2 . This probability may be regarded as giving a 
measure of the efficiency with which the system produces and detects a neutrino of a given 
mass. For convenience, the probabilities have been normalized to their values at m=0. In 



if iVe<Sl, that is, if the probability of detecting the neutrino in the bulk detector is much less than 
one (which is certainly the case). In the case at hand suppose that t\ corresponds to a time before 
any appreciable flux has arrived at the detector and ijv=ii + T to a time after all of the flux has 
passed. Then 

N (N — 1) N (N — l) f tl+T 

Y,V(x Dl ti)= y J - J2'P(xD,t i )At~ ± ' / dt D P(x D ,t D ). 

i=i i=\ Jtl 



11 



each case, the solid line gives the exact result and the dashed line shows the approximation 
for non-threshold masses derived in Eq. ( j2.34| ). Clearly the approximation is quite good if the 
mass is not too close to the neutrino production and detection thresholds. Furthermore, it is 
clear that this detector can be made "ideal" (that is, the probability to detect a neutrino may 
be made mass- independent) by using suitably small values for a Xl 2 . The mass-dependence 
for large a Xl 2 occurs for the same reason as in the case of the step function detector and is 
due to the fact that the source and detector in our model cannot recoil (see the discussion 
following Eq. ( |2.30| )). The dash-dotted line shows a plot of l/v(m) for comparison. This 
would be the analogous efficiency found in a wave packet calculation || . Since our detector 
may be made mass- independent by a suitable choice of the a Xi , we see that a "well designed" 
detector will not exhibit such an effect. 



B. Several Neutrinos 

Now that we have studied the characteristics of the source/detector system in the single- 
neutrino case, we turn to the case in which there are several neutrino fields coupled to the 
source and detector. Suppose that there are N different neutrino mass eigenstates. Then, 
in order to model the real-life situation, we suppose that there are also several different 
types of sources and detectors, each of which couple to a "weak eigenstate" which is a given 
unitary linear combination of the neutrino mass eigenstates. The action of Eq. ( [2.1| ) is then 
generalized to 

S = j d A x (£° + C- mt ) + J dtL° q , (2.36) 

where 

£$ = -£^)( D + ™?)<M x )' (2-37) 

i 

^ = E [#*(*)??(*) - 0?<tf(t)£(t) + q?(t)q«(t) - nfq?m(t)} , (2.38) 

a 

Ant = -E[ e iW (W^J(a09i(*)Mx) +^<M*)9i at (*K(x)) 

+e 2 (t) (W: i 0j(x)g 2 °(t)/ i2 (x) +W m i (x)g 2 Qt (t)^(x))] , (2.39) 

and in which U is a unitary matrix. Note that the subscripts "1" and "2" on the functions e 
and h and on the fields q refer, respectively, to the source and detector. These should not be 
confused with the subscripts on the fields <fi which refer to the mass eigenstates. Also note 
that we have taken e and h to be independent of the flavour or mass eigenstate in question. In 
principle there could be such a dependence, but including it would unnecessarily complicate 
our analysis. In what follows, we shall also set 0"=f2j, for all a, in order to "idealize" our 
sources and detectors. 

The experimental set-up which we wish to consider is a simple generalization of that 
given in the previous section. In this case we imagine that the initial state of the system 
has an a-flavour "source" oscillator in its first excited state and that the final state has a 
/3-flavour "detector" oscillator in its first excited state. The amplitude for this process may 
then be calculated as in the single-neutrino case and we find 
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= - £W/9iW^(0;0; lp\ J dt'dt"d 3 x'd 3 x"e 1 (t')e 2 (t") 

^,■(^(0^*)?? (*>iC«OI°; 1«; 0) 
= -J2 U PiKi[ dt'dt"d 3 x'd 3 x''dk i ei{t')e 2 {t'')h 1 (x')h;{x") 

i 

x exp - fi a )t" + i(E t - fix)*' + zk • (x" - x')] , 

in which we have defined 



Ei = Jk 2 + m 2 . 



[2A0) 



(2A1) 



[2A2) 



Taking hi, h 2 and e\ to be gaussians with widths a Xl , a X2 and as in the single- neutrino 
case (see Eqs. ( |2.18| ), ( |2.19| ) and ( 2.20 )), we may further simplify this expression 



A 



2ne\o tx 
4ti 2 xd 



Y,Uf*Ki dt"e 2 (t") dEexp[-z(E-n 2 )t" 

■ J —oo Jrrii 

sin(/cjX£)), (2.43) 



-liE-nifal-lk^l 



where 



E 2 — mf. 



(2.44) 



This expression is clearly just a simple generalization of Eq. ( 2.21| ). The final step in our 
calculation is to substitute the expressions fl2.23| ) and ( [2.24|) for e 2 (t) in the step function 
and gaussian detector cases. This yields 



C^fe-*i)E¥» dE 



°° sin [{E-n 2 ){t 2 -ti)/2\ 



sin(kiXn) 



x exp 



[(E-n 2 )(t 2 -ti)/2\ 

-~{E - n 2 )( tl + t 2 ) - \{E - ^fal - + al 2 



Af- = NV2^a t2 J2^K 



dEexp 



-~kHcr 2 Xl +(Tl 2 )-i(E-n 2 )t D 



X2' 



l - {E -^ifal-\{E-^ 2 fal 
sin(/ciX D ), 



(2.45) 



(2.46) 



where N is defined in ( j2.27| ). 

We are finally in a position to define the oscillation probability as a function of distance 
for the two cases. In both cases our definition of the probability is a "physical" one. We 
imagine that the source produces neutrinos of type a (a=e, //, r, . . .) and that we set a 
/3-neutrino detector at some distance xd from the source. We prepare the source (or an 
ensemble of identically prepared sources) in an excited state, wait a long period of time, and 
then check to see if the detector has been excited. After repeating this experiment enough 
times to get good statistics, we repeat the procedure with a /^'-neutrino detector, and so on. 
The probability to observe a (3 neutrino is then simply the number of events observed in 
"/9-mode" divided by the total number of events in all modes. Since we have attempted to 
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make our source/detector system as "ideal" as possible, there are no further corrections for 
detector efficiencies or other effects of that nature. The normalized coherent and incoherent 
oscillation probabilities may then be defined as 



A s ^Jx D , ti,t 2 ) 



(2.47) 



V, 




+ t 2 ' 

J °° dt D AfZ S p (xD,tD, eta ) 



(2.48) 



Y,/3' fo° dt D A z ^p,{x D ,t D ,a t2 ) 



2 



It is understood in the first expression that t\ is taken to be some time before the first bit 
of neutrino "flux" arrives at the detector. 

The expressions which we have derived for our two types of detectors are in forms which 
are amenable to numerical calculation. The coherent probability may be found after a single 
integration over energy and the incoherent probability requires two integrations, one over 
energy and one over time. In the two-neutrino case, the time integral in Eq. (|2.48| ) may be 
done by hand, but this is not possible in general for more neutrinos. The reason for this is 
that the integrand is no longer symmetric under £d— >— to due to the possible presence of 
phases in the mixing matrix U. 

Let us examine the case with two flavours in some detail. In this case, the matrix U may 
be taken to be a real orthogonal matrix parametrized by one angle, 9. The time integral in 
the numerator of ( |2.48| ) may be performed explicitly and we find 



Fig. H shows several plots of the flavour-conserving probability V e ^ e (x£)) as a function of 
xd for two relativistic neutrinos, using both the "coherent" and the "incoherent" detector. 
The various parameters chosen for the plot are as indicated in the figure. Recall that Qi and 

(set equal here) are the energies of the excited source and detector, respectively. Since 
we have chosen to set cx^-Ccr^, the energies of the mass eigenstates are approximately equal 
to Q and their momenta are determined by their energies. The values employed here for 9, 
mi and m,2 are chosen merely for the purpose of illustration. Note that the curves do not 
go all the way to Xd=0, since our formulas are not valid for very small Xp-f\ 

Figs. 0(a) and (b) show plots of the probability for detecting the same-flavour neutrino 
as was emitted in the case of an "incoherent" detector (see Eq. ( |2.48| )) for several different 
values of the time resolution of the detector, a t2 . The dotted curve in (b) is the analogous 



9 Recall that we require the source to turn off before the detector turns on in order that we may 
drop one of the time-orderings in the neutrino propagator. The reader is referred to the discussion 
following Eq. ( |2.15| ) for more details on this point. 




x exp [-(E - n x ) 2 a 2 h — (E — fi 2 ) 2 o£ 



(2.49) 
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result derived using the wave packet approach ||. This appears to be a good approximation 
to our result in the limit as o t2 ~ >0. It is clear from these plots that the coherence length 
of the oscillations is dependent on the time resolution of the detector; that is, as was noted 



in Ref. [fLql , a long coherent measurement in time is capable of "reviving" oscillations of 
neutrinos whose mass eigenstate wave packets have become physically separated. This effect 
is particularly striking in the case of the probability detected by the coherent detector, shown 
by the solid curve in Fig. 0(c). In this case the oscillations appear to have been completely 
revived even after, according to an "incoherent" measurement (dotted curve), the wave 
packets have completely separated. 

We have already discussed to some extent how it is possible for a long coherent mea- 
surement in time to revive the oscillations of neutrinos even after the mass eigenstates have 
separated spatially. Essentially, the accurate measurement of the energy picks out the plane 
wave in the wave packet which has existed coherently through both pulses Our present 



approach allows for a complementary way to view the situation. The question of whether 
the wave packets corresponding to two mass eigenstates have separated or not depends on 
the temporal and spatial resolution of the detector. We may demonstrate this effect by way 
of an example. Let us take a "snapshot" of the wave packets corresponding to two different 
mass eigenstates at a fixed time t£>=150 using the incoherent detector with different widths, 
a t2 . Figs. |3|(a) and (b) show the detection probabilities (given by Eq. ( |2.31|) , but separately 



normalized over xp) for the two mass eigenstates. In (a), the time resolution of the detector 
is taken to be cr <2 =l and there is almost no overlap between the two wave packets. Indeed, 
comparison with Fig. ||](a) shows that, for xd~150, the oscillations have been almost com- 
pletely damped out. If the detector is taken to have a broader time resolution as in Fig. ||(b), 
however, the wave packets appear to have a non-negligible overlap. In this case the width 
due to the finite time resolution of the detector has been added to the original widths of the 
wave packets. Comparison with Fig. 0(a) shows that in this case the oscillations have not yet 
been wiped out for xd~150. From this point of view, then, the fact that the conventional 
"wave packet" approach for relativistic neutrinos agrees with the source/detector approach 
(see Fig. ||](b)) for very small a t2 is not that surprising. The wave packet approach simply 
ignores the finite time resolution of the detector. 

The "coherent" and "incoherent" probabilities, Eqs. (|2.47| ) and (|2.48|) , may both be reli- 
ably approximated in the relativistic limit. Setting Q=Qi=Q<2 for convenience, we obtain 10 

V c a %(x D ) ~ -^YM^^M^i exp fe - kj)x D + (m 2 + m)){a 2 xi + oJ,)/2] (2.50) 



and 



Va™p{xD, o- t2 ) ~ ^-.y^U MM ijlLt exp [i(ki - kj)x D + (m 2 + m]){a 2 Xl + a 2 2 )/2 



10 We have used the approximate form of ,4 ste P given in Eq. ( p. 28 ) in order to derive Eq. ( p.50p . 
Also, in deriving Eq. ( 2.51| ) we have dropped the highly oscillatory terms in the integrand since they 
are strongly damped for xr>>&ti + <T *2- Recall that our calculation is only sensible for XD>&t 1 + &t 2 - 
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x exp 



X 2 D (1/Vj - 1/Vj) 



2l 



for the coherent and incoherent cases, respectively, where we have defined 



2/2 i 2 \ 



(2.51) 



(2.52) 



n 2 



ki 



Vi = ki/Vt. 



mi 



^2.53) 
(2.54) 



These expressions are identical except for the damping of the cross-terms which occurs in 
the approximation for the "incoherent" case, Eq. ( |2.51| ). Note that the oscillation length 
which may be extracted from either of these expressions is exactly what one finds in the usual 
approach, with no spurious factor of "2." The approximation for the "coherent" case contains 
no damping whatsoever, demonstrating that an infinitely long coherent measurement does 
indeed completely revive the oscillations of the neutrinos! We also note that, while our 
expression for the incoherent case, Eq. ( |2.51| ), bears some resemblance to the analogous 
expression obtained in the wave packet approach (see, for example, ||), our expression 
has an intrinsic dependence on the temporal and spatial resolution of the detector which 
is ignored in the wave packet approach. (This dependence on the detection process has 
also been investigated recently in P2"| .) Finally, we note the absence of factors of 1/vi 
pre-multiplying the exponentials such as can occur in wave packet calculations. 



C. The Non-relativistic Case 

It is worthwhile to consider briefly the oscillations of non-relativistic neutrinos in our 
toy model. Let us assume that one of the mass eigenstates is relatively light and let us 
study the behaviour of the oscillation probability as the mass of the other neutrino (in the 
two-neutrino case) is varied. Furthermore, let us restrict our attention to the case of the 
incoherent detector, which is the more realistic of the two detector types. As the mass of 
the heavier neutrino increases, the packets separate more quickly and, for sufficiently non- 
relativistic neutrinos, the oscillations are damped out almost immediately. It is convenient, 
then, to simply study the asymptotic expression 



^%K) = liva^V^ixu, a t2 ). (2.55) 



The main non-relativistic effect in our toy model is due to the model's dependence on a xi 2 
rather than being related to nonrelativistic effects of the oscillations themselves. Recall from 
our discussion in Sec. [11 A| that our source and detector are more efficient at producing and 



detecting non-relativistic neutrinos (see also Fig. |I[) This dependence skews the results for 
the oscillations, as one might expect. 

In Fig. |] we have plotted the probability for a v e to be detected as a v e in the limit 
as xd— >oc (V^ e in Eq. Q2.55p ) as a function of the mass of the heavier neutrino. The 
various curves correspond to different values of cr Xl 2 , the spatial resolution of the source and 
detector. For larger values of a Xl 2 , this probability is indeed skewed quite dramatically due 
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to the fact that the heavier mass eigenstate starts to dominate the probability distribution. 
Recall our earlier explanation as to why this occurs in our model. Since our source and 
detector do not "recoil" when the neutrino is emitted or absorbed, the upper limit on 
the neutrino's momentum is given by fc max ~l/cr xi a (the reader is referred to the discussion 
following Eq. ( |2.30| )). We emphasize, however, that this effect is an artifact of our model and 
would not be expected to occur in more realistic models. We shall discuss this point further 
below when we consider how our approach might be extended to the more realistic case in 
which the neutrino is not the only decay particle emitted. Note also that as m 2 increases 
above the production/detection threshold all of the solid curves approach the same value of 
cos 2 9. (How abrupt the threshold is depends on how large a tl and a t2 are, of course.) The 
dotted curve shows, for comparison, the result which is obtained if the detection efficiencies 
for the mass eigenstates are weighted by We see no evidence in our model for this type 
of behaviour. 



III. TOWARDS A MORE REALISTIC CALCULATION 

In this section we show how the bosonic model of the previous section may be modified 
to account correctly for the fermionic nature of the neutrinos (which we shall assume to 
be Dirac neutrinos) and for the V — A nature of neutrino interactions. This will have the 
added benefit of prefering neutrinos with nonzero momentum. Once again the source and 
detector will be modeled by harmonic oscillators. This time, however, the oscillators will 
be coupled to the usual V — A leptonic current rather than simply to the neutrino field. As 
a result, the interactions at the source and detection points will involve both the neutrino 
and its associated charged lepton. It is convenient to take the initial state to consist only 
of the source and detector, both in their first excited states. The source decays by emitting 
a neutrino and its associated charged anti-lepton, and the detector decays by absorbing the 
neutrino and emitting another charged lepton: 

(source)* — > v(k) + (pi) + (source) 

^ u(k) + (detector)* -> lp(p 2 ) + (detector). (3.1) 

This sequence of events is illustrated schematically in Fig. |^. The system may be described 
by the following action 

S = Jd'x (C° v + Ant) + / dtL° q , (3.2) 



where 



£° u = ^2vi(x) (i$ - rrii) Vi(x), (3.3) 

i 

L ° q = J2 [#*(*)#(*) - V 2 iq ?\t)q?(t) + $(t)q«(t) - tf 2 q? \t)q«(t)] , (3.4) 

a 

Ant = -£ [ei(t)qf (t)fci(x)j2(x) + e 2 (t)q«(t)h 2 (x)j°J(x) + h.c] , (3.5) 
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and where j^{x) is the zerothQ component of the leptonic V — A current 

f a {x) = J2K^YPlL(x), l a = e,fi,r,..., (3.6) 

i 

with Pl=(1— 7 5 )/2. Once again, eim and hm) are functions which parametrize the temporal 
and spatial couplings of the neutrino and lepton fields to the source (detector). 

The calculation of the amplitude proceeds in complete analogy with the calculation for 
the bosonic case and we shall omit most of the details. As above, we take ei(£) {^{t)) to 
be a gaussian of width a tl (&t 2 ) centered at t=0 (t=t£>) and hi(x) (/i 2 (x)) to be a gaussian 
of width a xi (a X2 ) centered at x=0 (x=x D ) (see Eqs. ( gTg ), (|CT| ), ( |CTD and (gg§). 
Thus, we omit here the case of the (coherent) "step function" detector and consider only the 
(incoherent) "gaussian" detector. Also, recall that the energies of the source and detector 
are fli and Q 2 , respectively. The amplitude to detect a neutrino of flavour (3 given that a 
neutrino of flavour a was emitted at the source is then given by 



r d?h r 1 

A a ^ = (27r)e?eK^ 2 £?VC J exp --(^ - E( Pl ) - E t fa\ 



d 3 k 

1 1 1 



(2tt) 3 2^ 

2„2 n. i „ |2_2 ■' ii „ 2 _2 



2 (fi 2 + ^ - S(p2))X - - |k + pxlX - - |k - p 2 | 2 < - iEit D + ik-* D 
x %(p2)7 ^(^ + ^)7°^(Pi), (3-7) 

in which the subscripts on the u and v spinors refer to their flavours; the spinors also have 
an implicit spin index which has been omitted. 

The above expression for the amplitude is qualitatively similar to the analogous expres- 
sion, Eq. (|2.46| ), derived previously in the bosonic model, with a few notable exceptions. On 



a technical note, we see first that it is no longer possible to perform the angular parts of the 
k integral exactly as was done in the previous case. This occurs because of the presence of 
the momenta of the charged leptons, pi and p 2 , which complicate the integrand somewhat. 
A related point is that now the neutrinos' momenta are not centered around zero, as was 
the case above. Rather, we have for the momenta 



and for the energies 



k^- Pl , (3.8) 
k « p 2 (3.9) 



Vti ps E(p\) + Ei, (3.10) 
Ei + ShxtEfa), (3.11) 



where Ei is the energy of the i th neutrino mass eigenstate. The relations ( |3.8|) -( [3~~i"l~l) are only 



approximate equalities since the degree to which each of them holds is determined by the 



11 In a more realistic calculation, one might perhaps couple the V — A current to a current repre- 
senting the initial and final nucleus . If these nuclei are sufficiently non-relativistic then it is a 
good approximation to consider only the zeroth component of the current. 
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relative sizes of a Xl , . . . , <Jt 2 • The fact that the neutrinos' momenta are not centered about 
the origin is rather encouraging because it indicates that this model would not be expected 
to have the (unphysical) feature that it favours non-relativistic neutrinos, as was the case 
in the bosonic model of the previous section. The final difference, compared to the bosonic 
case, is the presence of the matrix element, up . . .v a , which contains all of the information 
regarding the neutrinos' spins. It is interesting to note the presence of the factor 

(3.12) 

which arises in this case in part due to the sum over spins of the neutrino u spinors, 



^2 s u s (ki)u s {ki). This same factor appears in the field theoretic calculation of Ref. [19 
but in that case is due to an integral in the complex fc plane which extracts the pole of the 
propagator |24| . We need not do any such integration since we always insist that our source 



be turned "off" before our detector is turned "on." This forces the neutrinos to always be 
on-shell. 

It would be possible at this point to proceed as we did in the previous section. First we 
could examine the response of the detector to the source by looking very carefully at the 
case in which there is only one neutrino. Armed with this knowledge we could define the 
probability in analogy with the bosonic case and study its behaviour as a function of the 
various parameters of the theory. While this progam might be deserving of future study, for 
now we shall content ourselves with a more qualitative examination of the generic features 
of this model. 

As we have noted, two of the main qualitative differences between this model and our 
former bosonic model are the different energy-momentum conservation equations and the 
presence of the matrix element in the integrand. A further difference is that in order to 
obtain the oscillation probability, we now need to integrate over the momenta of the two 
outgoing charged leptons. Since the k integral in the expression for the amplitude is expected 
to be dominated by values of k which are parallel to x B f20fl , the pi and p 2 integrals would 
similarly be dominated by values anti-parallel and parallel, respectively, to x^, due to the 
damping terms in the exponential of Eq. (|3.7| ). In order to get some idea of the effect of the 
matrix element as a function of the neutrino's mass, then, let us evaluate it when all of the 
momenta are parallel (or anti-parallel) to xq. (A somewhat similar analysis to the following 
may be found in Ref. p5| .) Choosing an explicit representation for the gamma matrices and 
adopting the normalization conditions of Itzykson and Zuber [F2(| pp. 57, 145-6, 201], we 



find that only two of the four helicity combinations of the leptons survive, yielding 

M+Umii = -IE, - k) <f (Pl) + '"° ti^EM±^i -»») (3.13) 

2 [Am a mp(E( Pl ) + m a )(E(p 2 ) + mp)} L/2 

ka— ( \ (t? i u\ (E( Pl ) + m a - Pl )(E(p 2 ) + 171(3 +p 2 ) 

M a ^p(mi) = -{Ei + k)— — — — — -Y72, 3.14) 

where fc=|k|, etc., and where the "++" and " " superscripts refer to the helicities of the 

lepton and anti-lepton. In the limit as the neutrino mass goes to zero, only the combination 
in which both leptons have negative helicity survives, since the exchanged neutrino can only 
have negative helicity in that limit. For non-zero masses it becomes possible to also produce 
lepton pairs with positive helicity. 
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The quantities which will occur in the oscillation probability are the squares of the matrix 
elements. Let us define 



hUpfa) = \Mitp(m)\ 2 /\M~e(0)\ 2 , (3.15) 
h-^(m) = \M-- p (m t )\ 2 /\M-- p (0)\ 2 . (3.16) 

Then h + (h~) gives some measure of the probability that the source/detector interaction 
gives rise to two leptons with positive (negative) helicity. Since the efficiency of the system 
at producing and detecting neutrinos of a given mass is determined to some extent by the 
functions h , it is useful to plot them as a function of the mass of the exchanged neutrino. 
It turns out that the energy-momentum conservation equations, Eqs. ( |3.8| )-( ]3~TT| ), are 



over-complete. Thus, for given values of the charged lepton and neutrino masses, for exam- 
ple, fli and Q2 may be found such that all of the conditions are met, but when the neutrino 
mass is varied, at least one of the conditions needs to be violated. This problem is related 
to the difficulty which occurred in the bosonic model (where momenta close to zero were 
favoured) and has its root in the fact that our source and detector are fixed and do not 
recoil. For the purposes of our plot, let us require that Eqs. (|3.8|) , (|3.10|) and (|3.11|) hold 



exactly - so that energy and momentum are conserved at the source and energy is conserved 
at the detector - and allow the momentum conservation at the detector, Eq. (|3.9| ), to be 
violated. As in our previous model, this can again be allowed by setting a X2 to be somewhat 
small.Q For the plot let us take a=(3= e, so that both the source and detector are sensitive 
to electron neutrinos. We then set 

Q 1 = 0.6 MeV, Q 2 = 0.5 MeV, 

m a = m p = me = 0.511 MeV. (3.17) 

Fig. H shows a plot of h^ e {m) and h~^ e {m) as a function of the neutrino mass. The 
"threshold" in this case is determined by the condition Qi=m e + m, where m is the neutrino 
mass. The upper curve corresponds to the negative helicity case and approaches unity as 
m— >0. The lower curve disappears in the same limit. For neutrino masses closer to threshold, 
fairly substantial deviations from the m=0 case are observed to occur. 

The plot in Fig. ^| should of course be treated with some caution, since it shows only the 
square of the matrix element evaluated at some "optimal" energy and momentum configu- 
ration. In general, the oscillation probability will also receive contributions due to energy 
and momentum configurations which are non-optimal. Furthermore, it has been found that 
the procedure which we have followed can lead to non-sensical results if the neutrino mass 
is taken to be large compared to the lepton mass.|^] In any case, however, the plot does 



On physical grounds we would prefer to allow momentum conservation to be violated somewhat 
rather than energy conservation. The reason for this is that in the former case, the small value 
required for a x is still of a reasonable magnitude compared to nuclear scales (it is on the order of 
several hundred fm in the example considered here), but the value which would be required for at 
would be far too small compared to any time scales in the physical problem. 

13 This occurs because, in our prescription, k and p 2 need not be the same. For very heavy 
neutrinos this starts to cause problems in this approach. 
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demonstrate something which might be regarded as "typical" : for non-relativistic neutrinos 
there will be a non-zero probability to produce charged leptons in the final state which have 
the "wrong" helicity configurations. Thus, particularly if the spin of the leptons were to be 
measured in a certain experiment, one could expect there to be quite strong mass effects for 
non-relativistic neutrinos. In our case, for example, there is a suppression of the negative 
helicity final states for large mass and a mild enhancement of the positive helicity ones. 

Since in this model the neutrinos no longer have their momenta centered about the 
(unphysical) value of "zero," one would expect in this case that the non-relativistic neutrinos 
would not be favoured, as was found to be the case in the bosonic model studied above. In 
fact, it is possible that there would be a suppression for non-relativistic neutrinos due to the 
phase space suppression of the final state leptons, for small momenta. This question could 
really only be answered by performing a thorough numerical analysis of the model, which 
we shall not do at this time. 

IV. DISCUSSION AND CONCLUSIONS 

Most phenomenological work in the field of Particle (z/, K, B, ...) oscillations describes 
the oscillations as a function of time and then converts the time dependence of the results to 
a space dependence. There have been many attempts in the literature to improve on these 
calculations by explicitly including the spatial dependence of the wave function. These 
approaches have necessarily lead to the description of the wave function as a wave packet. 
It has been shown that several recent claims that such wave packet approaches lead to 
different results than the simple time-oscillation approach, are incorrect and that a proper 
wave packet calculation leads to the "expected" results. 

In this paper we have presented a novel approach to the study of the spatial dependence 
of neutrino (and other particle) oscillations. We have done this by coupling the neutrino field 
to an idealized, localized model of a source and detector which we have chosen to describe as 
simple harmonic oscillators which can be excited or de-excited by the absorption or emission 
of a neutrino. The system begins with the source in the first excited state and the detector 
in its ground state. We then compute the probability that at a much later time the source is 
in its ground state (so that it has emitted a neutrino) and the detector is in its first excited 
state (so that it has absorbed a neutrino). This probability is evaluated as a function of 
the distance between the source and the detector and it depends, in detail, on the spatial 
extent of the source and the detector as well as on the length of time for which each is on. 
We have seen how to use this dependence to obtain a better understanding of how neutrino 
oscillations depend on the time resolution and the coherence properties of the source and 
the detector. We have also seen how our approach is useful in clarifying several subtle issues 
related to the Quantum Mechanics of neutrino oscillations. 
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APPENDIX A: APPROXIMATE AMPLITUDE FOR THE COHERENT 

DETECTOR 



In this appendix we shall derive an approximation for the t 2 -^oo limit of the integral 
given in Eq. ( J2.25| ) and investigate under what circumstances the approximation is valid. 

The form for the integral given in Eq. (|2.25 ) is convenient for numerical work, but is not 
particularly convenient for the limit which we wish to consider. Let us instead go back to 
the definition of this expression, gotten by inserting Eq. (|2.23|) into Eq. ( p. 21 ). We may now 
formally take the limit as t 2 ^oo by giving Q 2 a small imaginary piece. This yields 
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where the limit e^0 + is understood. This integral may be simplified by employing the 
relation 
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to obtain 
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In order to approximate Eq. ( |A3|) it is useful to make a change of variables. On the interval 
(m, Q 2 ) we define E=Q 2 — E and on (Q 2 ,oo) we define E=E — Q 2 . Then the integral in 
may be approximated by 
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where Af2=f2 2 — ^1 and where the only approximation so far is that the interval (£7 2 , oo) has 
been truncated to (f2 2 ,2f2 2 — m). This approximation is valid if the major contribution to 
the integral comes from energies close to f2 2 . in order to further approximate the integral, let 



us make the ansatz that the integral in (|A5|) is dominated by values so close to E=0 that is 
valid to set E=0 in the gaussian pieces. At the end of the calculation we will be able to see in 
which cases this is a reasonable approximation. When dealing with the oscillating terms we 
must be a bit more careful. Writing the sine's in terms of exponentials and Taylor-expanding 
the arguments to first order in E (which essentially amounts to ignoring the spreading of 
the wave packets) leads to the following approximation for (|A5|) 
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where 
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The final step in the approximation is to note that, if 



— ± h > a h 
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(A8) 



and if |Afkr tl | is less than or of order unity, then we may approximate the sine terms by 
delta functions, since 



L^oo x 

This brings us to the desired result 
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Note that the condition in Eq. ( |A8| ) simply requires that the detector be turned on before 
any appreciable amount of flux reaches it. 
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FIG. 1. Plot of the "incoherent probability" "Pi nco h (Eq. ( |2.33| ), normalized to its value at m=0) 
as a function of mass for the case of a single neutrino, taking a Xl 2 =0.1, 0.05. In each case, the solid 
line shows the exact result and the dashed line shows the result obtained in the approximation of 
Eq. (|2.34| ). The dash-dotted line shows a plot of \/v(m) for comparison. 
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FIG. 2. Oscillation probabilities as a function of distance. The two curves in (a) correspond 
to the "incoherent" detector with time resolutions o"t 2 =l, 2. The solid curve in (b) gives the 
"incoherent" probability for cr t2 =0.1. The dotted curve shows the analogous result obtained in the 
wave packet approach. The solid curve in (c) shows the probability measured by the "coherent" 
detector. 
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FIG. 3. "Snapshots" of two mass eigenstate wave packets using incoherent detectors with dif- 
ferent time resolutions. The wave packets have been individually normalized over X£>. In (a) the 
time resolution of the detector is such that the wave packets appear to be nearly separated, while 
in (b) the same wave packets appear to overlap due to the broader temporal resolution in that 
case. 
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FIG. 4. Plot of the constant flavour-conserving probability, V^ e , defined in Eq. ( [2.55 ) as 
a function of the mass of the heavier neutrino. The solid lines correspond to spatial widths 
a xi 2 =0.01, 0.05, 0.1, 0.15 and the dotted line shows the value obtained in the wave packet approach 
if the contributions corresponding to the various mass eigenstates are weighted by 1/vi. 
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source detector 




FIG. 5. A schematic illustration of the sequence of events in the source/detector system for 
fermionic neutrinos considered in Sec. [TO] . The excited source decays by emitting a neutrino and 
its associated anti-lepton. The detector subsequently absorbs the neutrino and emits a lepton. 



29 



1 rr 



0.8 



0.6 



0.4 



0.2 







n, = 0.6 MeV 



Q 3 = 0.5 MeV 



rn = 0.511 MeV 



J I L 



I I I 







20 40 60 

Neutrino Mass (keV) 



h- 



80 



FIG. 6. Plot of the two functions hf_^ e and h~_^ e as a function of the neutrino mass. These 
provide a measure of the probability to produce lepton pairs with helicity +1 and —1, for h + and 
h~, respectively. 
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